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Abstract 

We present two iterative methods for computing the global and partial 
extinction probability vectors for Galton- Watson processes with countably 
infinitely many types. The probabilistic interpretation of these methods 
involves truncated Galton- Watson processes with finite sets of types and 
modified progeny generating functions. In addition, we discuss the connection 
of the convergence norm of the mean progeny matrix with extinction criteria. 
Finally, we give a sufficient condition for a population to become extinct 
almost surely even though its population size explodes on the average, which is 
impossible in a branching process with finitely many types. We conclude with 
some numerical illustrations for our algorithmic methods. 

Keywords: multitype branching process; extinction probability; extinction 
criteria; iterative methods 

2010 Mathematics Subject Classification: Primary 60J80 

Secondary 60J05;60J22;65H10 



* Postal address: Department of Mathematics and Statistics, The University of Melbourne, VIC 3010, 
Australia, sophiemh@unimelb.edu.au. 

** Postal address: Departement d'Informatique, Universite libre de Bruxelles 1050 Brussels, Belgium. 
latouche@ulb.ac.be. 

*** Postal address: School of Mathematical Sciences, The University of Adelaide, 5005, Australia, 
giang . nguyen® adelaide . edu.au. 



1 



2 



S. Hautphenne, G. Latouche, and G. T. Nguyen 
1. Introduction 



Branching processes are powerful mathematical tools frequently used to study the 
evolution of collections of individuals over time. In particular, multi-type Galton- 
Watson processes represent populations in which individuals are classified into different 
categories and live for one unit of time. Each individual may reproduce at the end of 
its lifetime, with reproduction rules dependent on its type. 

When the number of types is finite, one extinction criterion is based on the spectral 
radius sp(M) of the mean progeny matrix M, the elements My of which are the 
expected number of direct offsprings with type j for a parent of type i. Moreover, the 
extinction probability vector q is the minimal nonnegative solution of the fixed-point 
equation q = P(q), where each component qi is the extinction probability given the 
initial type i, and P(-) is the progeny generating function of the process. Harris [3] and 
references therein present a comprehensive analysis of extinction criteria and extinction 
probability for Galton- Watson processes with finitely many types. 

To allow, as we do here, the set of types to be infinite gives rise to three main 
challenges. First, as the mean progeny matrix M has infinite dimension, one has to look 
for a replacement to the spectral radius as an extinction criterion. Second, one needs to 
determine how to compute the extinction probability vector q which now has infinitely 
many entries. Third, the concept of extinction has to be defined carefully: when there 
are infinitely many types, it is possible for every type to eventually disappear while the 
whole population itself explodes. We use the term global extinction to indicate that the 
whole population becomes extinct, and represent by q the probability vector for this 
event; we refer to the event that every type becomes extinct as partial extinction, and 
denote its probability vector by q, with q < q, naturally, and the question is whether 
they are equal or not. 

Galton- Watson processes with infinitely many types have been much investigated 
already. Moyal [10J assumes that the types belong to an abstract space and proves that 
the extinction probability is a solution of the fixed point equation s = P(s). Mode 
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Theorem 7.2], for a restricted family of progeny densities, gives an extinction criterion 
based on the spectral radius of some integral operator. Focusing on denumerably 
infinite sets of types, Moy O |H| and Spataru [T3] use ergodic properties for infinite 
matrices, and analyse in special cases the role of the convergence norm of M as an 
extinction criterion. Recently, some authors in the literature of branching random 
walks have defined local survival, meaning that for every given type i and arbitrarily 
large epoch T there is at least one individual of type i alive at some time t > T , with 
global survival meaning that at least one individual is alive at any time, and strong local 
survival, when the two have the same probability. We refer to Bertacchi and Zucca [I], 
Zucca [16], and to Gantert et al. [2]. 

There is, however, no simple general extinction criteria for Galton- Watson processes 
with countably infinitely many types so far, and the question of actually computing 
the extinction probability vector has received scant attention, if any. 

Our main result is the development of two algorithmic methods for computing the 
global and the local extinction probability vectors q and q. The methods, which are 
presented in Section [31 have a physical interpretation based on two truncated Galton- 
Watson processes with finite sets of types. They may be applied to both irreducible 
and reducible branching processes with countably infinitely many types. 

In Section@]we discuss some extinction criteria expressed in terms of the convergence 
norm of the mean progeny matrix M in the irreducible case, or of irreducible sub- 
matrices of M when M is reducible. We also give a sufficient condition under which 
the population becomes extinct almost surely while its expected size tends to infinity. 
That condition implies that the asymptotic growth rate of the process may depend on 
the distribution of the initial individual's type. 

In Section [SJ we provide some numerical illustrations. Our examples are taken from 
two classes of processes for which the matrix M is tridiagonal (and irreducible) or 
super-diagonal (and reducible). 
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2. Preliminaries 



Consider the process {Z n = {Z n i,Z n2 , . . .)}„ £ b, where Z n e is the number of indi- 
viduals of type I alive at the nth generation, for t in the countably infinite set of types 
S = {1,2, 3,...}. Unless otherwise stated, the process starts in generation with one 
individual. 

We denote by Pij for j = (j\ , j 2 , . . .) the probability that an individual of type i gives 
birth to ji children of type 1, j 2 children of type 2, etc., and the progeny generating 
function P%{s) of type i e S is given by 

oo 

jeK°° k=l 

with s = (si,s 2 ,...), s t e [0,1] for all i. We define P(s) = (Pi(s), P 2 (s), . . .). The 
mean progeny matrix M is defined by 



for i,j e 5, 

s=l 



and Mij is the expected number of direct offspring of type j born to a parent of type i. 
The process {Z n } is said to be irreducible if M is irreducible, and it is reducible 
otherwise. 

The total population size at the nth generation is \Z n \ — Y^TLi Znl, an d we denote 
by ipo the type of the first individual in generation 0. The conditional global extinction 
probability vector, given the initial type, is q = (qi, q 2 , . . .) where 

qi = P[ lim \Z n \ = 0| <fo = i] for i e S. 

n— >oo 

This is the usual conditional probability that the whole population eventually becomes 
extinct, given the type of the initial individual, and we write that q = P[lim„_ i . 00 \Z n \ = 
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for short. The vector q is the minimal nonnegative solution of the fixed-point equation 



P(s) = s. (1) 

This equation has at most two distinct solutions, 1 and q < 1, if M is irreducible, and 
potentially infinitely many solutions otherwise (Moyal [10!, Spataru [13J). 

The conditional partial extinction probability, given the initial type, is q = (q%, q2, . . .) 
where 

<ji = P[W : lim Z n i — 0| ipa = i] for i 6 S. 

n— >oo 

In the irreducible case, Zucca [TB] observes that lim JWOO Z n t — for all types I if and 
only if the limit is zero for at least one type, regardless of the initial type. 

The vector q is also a solution of ([T]). Indeed, by conditioning on the progeny of the 
initial individual and using the independence between individuals, we readily obtain, 
for any i, 

qi = P[W : lim Z n£ = 0| cp = i] 

n— ►oo 

oo 

= Y PuTT p [ W: lim Z rU = 0\vo=k] jk 
= Pi(q)- 

When the set of types is finite, global and partial extinction are equivalent, but this 
is not necessarily the case when the set of types is infinite: by Fatou's Lemma: 



lim \Z n \ = lim Z n£ > lim Z n£ , 

n— ¥oo n— >oo * — • L — » n— voo 



so that 



P[ lim \Z n \ = 0|v? = i] < P[W : lim Z nl = 0|^ 
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for i,£ G S and < q < q < 1. 

As the vectors q, q and 1 are all solutions of (fTJ) and since there are at most two 
distinct solutions in the irreducible case, the following lemma is immediate. 

Lemma 2.1. If M is irreducible, and q < 1, then q = q. □ 

In the irreducible case, q = q < 1 is equivalent to strong local survival in the 
terminology of branching random walks and, although it is expressed differently, this 
sufficient condition is observed in [5] with the assumption that M is tridiagonal, and 
in [TB] for the general case. When M is reducible, it is possible that q < q < 1, and 
we give an example at the end of Section [5] 



3. Computational aspects 

In this section, we develop iterative methods to compute the extinction probability 
vectors q and q. The procedures apply for both irreducible and reducible processes. 
The underlying idea is to compute approximations of the infinite vectors q and q by 
solving finite systems of equations in such a way that the successive approximations 
have probabilistic interpretations: for q we use a time-truncation argument, and for q 
a space-truncation argument. 

3.1. Global extinction probability 

Denote by N e the generation number when the process becomes extinct. Clearly, 
q = ¥[N e < oo \ ipo\- Let Tk = {k + 1, k + 2, . . .} be the set of types strictly greater than 
k, and define the first passage time = inf{n : X^eT fc Z n £ > 0}, for k > 0. This is 
the first generation when an individual of any type in Tk is born. Furthermore, define 

q\ k) =V[N e <T k \ip Q = i\, 

the conditional probability that the process eventually becomes extinct before the 
birth of any individual of a type in given that the initial individual has type i, and 
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Lemma 3.1. The sequence {q^} k >a is monotonically increasing and converges point- 
wise to the global extinction probability vector q. 

Proof. Clearly, T k D T k+1 for all k, and r fe < T k+1 , so that [N e < T k ] C [N e < T k+1 ], 
and q( k > < q^ k+1 \ Therefore, for any i, 

lim q\ k) = lim P[7V e < r fe | <p = i] 

k— >oo k^-oo 

= F[N e < lim r k | (f = i] 

k—too 

= ¥[N e < oo | ip = i] 
= Qi, 

which completes the proof. 

(k) 

By definition, q\ ' = for all i £ T k . Consequently, 
q (0) = (0,0,...), 

qW = (q[ k \...,qi k) ,0,0,...) for k > 1. 

Thus, at the fcth iteration one only needs to compute the finite vector w/ fe ' = (q[ k , . . . , ), 
which we do as follows. Consider a branching process {W^}nEK which evolves like 
{Z n } under taboo of the types in T k . The taboo progeny distribution associated 
with types i e {1, . . . , k} in {wi^} is defined as 



h(h,...,j k ) — Pi(h,-,jk,0,0,...)- 



If the process is irreducible, then J2j^n k fi(j ± j k ) — 1' f° r 1 — * — ^' w ^ n at l eas t 
one strict inequality, and we need to add an absorbing state A to the state space 
M fe of {Wi k) } to account for the missing probability mass. Obviously, absorption in A 
precludes extinction, and is the vector of probability that {W^} becomes extinct 
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before being absorbed in A, given the type of the initial particle. In consequence, 
is the minimal nonnegative solution of the finite system of equations 

Si =F j; (fc) ( Sl , S2 ,..., Sfc ), forl<i<fc, (2) 

where F^ k \s) = Pi(s\, . . . , Sfc, 0, . . .) is the probability generating function of the 
distribution . 

We may compute u/ fe ) by linear functional iteration, for instance, on the fixed- 
point equation @: one easily verifies that, for any k > 1, the sequence {i«( fc >™) = 
(w[ k ' n \ . . . , toi* )} n >o recursively defined as 

(k,n) TTi(fe) / (k.n— I) (k,n— 1)\ r 1 ^ • ^ 7 

w 8 = F> (wj , . . . , u>£ ) for 1 < i < k. 

for n > 1, with u/ fc '°) = (0, . . . , 0), satisfies 

^(fc.n) = p[jv e < Tfe and AT e < n | <p ], 

and is, therefore, monotonically increasing to «/ fe ). In practice, we would terminate 
the functional iteration for a given k when becomes sufficiently 

small. 

3.2. Partial extinction probability 

Here, we associate to the branching process {Z n } a family of processes {Z^ — 
(Z^i , ■> ■ ■ OlneiNi f° r k > 0, obtained as follows: for a given k, we count neither the 
individuals of types in X^, nor any of their descendants, whatever their types. It is as 
if all individuals of types in 7\ became sterile. Define q^ k ' to be the global extinction 
probability vector of the process {Zn}. 

Lemma 3.2. The sequence of vectors {q^}k>o * s monotonically decreasing and con- 
verges pointwise to the partial extinction probability vector q. 
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Proof. Obviously, 

z^ = (zi%z^,...,zHo,o,o,...) 

W 7 (fc+1) 7 (fe+l) y(k+l) 7 (fe+l) nn \ 

^ l Z nl > Z n2 )---> Z nfc ' Z n(fe+1) ' U ' U ' • " ^ a ' S - 

so that, for n fixed and k —> oo, monotonically converges to Z n . Furthermore, 

lim |Z«|< lim |Z( fe+1 )|, 

n— >-oo n— »oo 

so that > q (fe+1) , for jfc > 0. 



Let = [linin^oo = 0] be the event that type I of {Z^} eventually becomes 
extinct, and let = C\i>\B^ be the event that all types of {zi^} eventually become 
extinct. We have 

gW=P[lim \Z^\=0\ Vo ]=V[A^\ V o], 

n—>oo 

since \zffi\ = ^iZ^\ contains only finitely many nonzero terms. Furthermore, 

B (fc+i) Q B (k) ) so that A(fc+ i) g i4(fc)j and 

A {k ) = f|[Iim Z„< = 0]. 

£>1 



Therefore, 



lim q (fc) = PfAco | po] = P[W : lim Z nt = | po] = 9, 



which completes the proof. 
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By definition of {Z^}, q^ k ' — 1 for all i £ Tk, and so 

g(°) = (l,l,...) ) 

g« = (#\.. .,^,1,1,...). 

To compute the finite vector — (q[ k ^ , . . . , % ), we may interpret {Zn^} as a 

Galton- Watson process with finitely many types and progeny distribution f\ defined 
as follows: 

jfe+i.ife+2,--->o 

and vft^ is the minimal nonnegative solution of the finite system of equations 

Si = F i (fc) (si,s 2 , • ■ • ,s k ), for 1 < i < k, 
where F^ k \s) — Pi(s%, . . . , s^, 1, 1, . . .) is the probability generating function of the 

~(k) 

distribution /y . That equation may be solved by functional iteration, as explained 
at the end of Subsection 13.11 

4. Extinction criteria 

When the number of types is finite and M is irreducible, it is well-known that 

• q < 1 if sp(M) > 1, 

• q = 1 if sp(M) < 1. 
If M is reducible, then 

• q ^ 1 if and only if sp(Af) > 1, 

where we write v ^ 1 to indicate that Vi < 1 for all i, with at least one strict inequality. 
Indeed, if M is reducible, there may exist some type (but not all) from which partial 
extinction is almost sure even if sp(M) > 1 (Hautphenne [4]). 
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We expect that, in the infinite countable case, some analogue of sp(M) also plays 
a role in determining if extinction occurs almost surely or not. This is the case for 
partial extinction, but not necessarily for global extinction. 

4.1. Partial extinction — M irreducible 

We denote by the mean progeny matrix of the process {Zi? ] } defined m 

Section [3721 and by the k x k north-west truncations of the infinite matrix M. 

As we do not count individuals with types in Tf., is given by 



M« = 















and it is clear that q {k) = (q[ k \ . . . , q[ k \ 1, 1, . . .) = 1 if sp(AfW) < 1, otherwise 
q {k) < 1, with (q[ k \ . . . , q (k) ) < 1 if is irreducible. 

We assume in this subsection that M is irreducible. The convergence norm of M is 
defined as follows. Let R be the convergence radius of the power series ^2 k>0 r k (M k )ij, 
which does not depend on i and j. The convergence norm v of M is 

v = R- 1 = lim {{MX} 1 /*; 

k— >oo 

it is also the smallest value such that there exists a noimegative vector x satisfying 
xM < vx (Seneta |12[ Definition 6.3]). Note that the convergence norm of a finite 
matrix is equal to its spectral radius. 

If we assume that all but at most a finite number of truncations M^ 1 are irreducible, 
then by Seneta ( |12[ Theorem 6.8]) the sequence {sp(Af( fe ))} is non-decreasing and 
linifc^oo sp(M( fc )) = and one immediately shows the following property. 

Proposition 4.1. Assume that M is irreducible. The partial extinction probability 
vector q is such that q < 1 if and only if v > 1. 

Proof. If v > 1, as sp(M^ fe ^) /• v, there exists some k such that sp(A/W) > 1 
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and such that is irreducible. Thus, (q[ k \ . . . ,<j^) < 1 and, since q^ \ q by 

Proposition 13.21 q < 1 in the limit. 

If v < 1, then for all fc, sp(A/W) < 1 and q (fc) = 1, which implies that q = 1. 

This is also observed in Gantert et al. [2] and Miiller [11] . who use different 
arguments. 

4.2. Partial extinction — M reducible 

Let us assume now that the matrix M is reducible. The sequence {sp(M^)} 
is still non-decreasing, but its limit might not be the convergence norm of M. Let 
v G [0,oo] denote the limit. The proof of the proposition below is very similar to that 
of Proposition 14. II and is omitted. 

Proposition 4.2. The partial extinction probability vector q is such that q = 1 if and 
only if v < 1, otherwise q ^ 1. 

In other words, there exists at least one type i such that q~i < 1 if and only if v > 1. 
The next question is, for which i does the inequality qi < 1 hold? We give below a 
necessary and sufficient condition for qi to be strictly less than one. 

We write i — » j when type i has a positive probability to generate an individual of 
type j in a subsequent generation, that is, if there exists n > 1 such that (M n )ij > 0. 
We define equivalent classes C\, C2, . . . such that, for each fc, if i G Ck, then j 6 CV- if 
and only if i — > j and j i, for all j. This induces a partition of the set of types S 
and we write that Ck ~ > Ci when there exist i € Cfc and j G CV such that i — > j. We 
denote by Mk the irreducible mean progeny matrix restricted to types in Ck, that is, 
Mk = (Mij)ij£c k j an d by Vk the convergence norm of Mk- 

Proposition 4.3. // i is a type in Ck, then the partial extinction probability qi is 
strictly less than 1 if and only if Vk > 1 or there exists a class Ce such that Ck — > Ci 
and vi > 1. 

Proof. Let i G Ck and assume that > 1. Then, by Proposition l4.1[ the probability 



Branching processes with infinitely many types 13 

that every type in Ck eventually becomes extinct, given that the initial type is in Ck, 
is strictly less than one, hence qt < 1. 

Now assume that i € Ck, and that there exists a class Cg such that Ck — > Cg and 
i>i > 1. Then, i — > j for all j £ Cg, that is, there is a positive probability that type 
i has type j £ Cg among its descendants; moreover, since v/_ > 1, starting from any 
j £ Cg, the probability that every type in Cg eventually becomes extinct is strictly less 
than one by Proposition 14.11 We thus obtain q\ < 1. 

If Vk < 1 and there is no class Cg such that Ck — > Cg and vg > 1, then all classes 
such that Ck — > Cg satisfy vg < 1. In other words, by Proposition 14.11 all the 
descendants of type i will generate a process which partially becomes extinct with 
probability 1. So partial extinction is almost sure if the process is initiated by type i, 
and we have = 1. 

4.3. Global extinction — M irreducible 

We assume again that M is irreducible. By Lemma 12.11 and Proposition 14. li we 
know that if v > 1, then q = q < 1, and if v < 1, then q < q = 1. One question which 
remains is to determine additional conditions that guarantee q = q = 1. 

The most precise answers are conditioned on the dichotomy property, which states 
that with probability 1 the population either becomes extinct or drifts to infinity 
(Harris |3J). In the finite case, this follows under very mild conditions but it is more 
problematic if the number of types is infinite. In particular, Tetzlaff [141 Condition 2.1 
and Proof of Proposition 2.2] gives the following sufficient condition for the dichotomy 
property to hold: it suffices that for all fc > 1, there exists an index rrik and a real 
number dj- > such that 

inf P[|^ m J = | tp =i, 1 < \Zx\ <k]> d k . (3) 

i 

This indicates that there is a positive, and bounded away from zero, probability for 
the population to become extinct. The next property is proved in |14| . 
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Proposition 4.4. Assume that the dichotomy property holds. If the limit liminf r n. 00 (Af 
is finite, then q = 1 . 

A direct consequence, which brings the convergence norm back into the picture, is the 
following. 

Proposition 4.5. Assume that the dichotomy property holds. If there exist A < 1 and 
x > such that xl < oo and xM < Xx, then q= 1. 

Proof. Under the assumptions of the proposition, xM n l < X n xl, which implies 
that lim n _>oo xM n l < oo. Applying Fatou's Lemma, we obtain xlim n _>oo(Af™l) < oo, 
which leads to \\m n ^>. 00 (M n l) < oo since x > 0. Thus, by Proposition 14.41 the result 
follows. 

If such a A exists, it is at least equal to v, and we remember that v < 1 is a necessary 
condition for q = 1. The difference between A and v, and the additional constraint 
imposed by this proposition, is that the measure associated to A must be convergent, 
which is not necessarily the case with the measure associated with v. 

4.4. Growth rate and extinction 

When the number of types is finite and the process is irreducible, the expected total 
population size increases, or decreases, asymptotically geometrically: independently of 
the initial type, E[|Z„|] ~ p n where p is the spectral radius of M . This is no longer the 
case when the number of types is infinite, and the evolution of E[|2„|] may depend on 
the distribution of (po. Actually, it is possible for a process to become globally extinct 
almost surely while the expected population size increases without bounds. This we 
show below, and we give one example in the next section. 

Assume that there exists a probability measure a.\ such that ol\M < Aicti with 
Ai < 1, and a probability measure oli such that ol^M > A2CI2 with A2 > 1- If, in 
addition, the dichotomy property holds, then q = 1 by Proposition 14.51 
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If (fQ has distribution c*i, then 
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E[\Z n \] = ai M n l<\? 

so that lim„_ ! . 00 E[|Z„|] = 0. In the contrary, if ipo has the distribution a.2, then by 
a similar argument limn-^ E[|Z„|] = oo and the extinction probability is equal to 
a 2 q = 1. 



5. Illustration 

We illustrate the results of the previous sections with two examples, one for which 
M is tridiagonal (and the process is irreducible) and one for which M is super-diagonal 
(and the process is reducible). 

5.1. Irreducible tridiagonal case 

This example corresponds to a homogeneous branching random walk on positive 
integers with a reflecting wall at z = 1. The mean progeny matrix is 



M = 



b c 

a b c 

a b c 



(4) 



where a and c are strictly positive, and b is nonnegative. 



Proposition 5.1. Assume that M is as given in Then, its convergence norm is 
v = b + 2^/ac, and there exists x > such that xM = Xx for all X > v. In addition, 
xl < oo if and only if X £ [v, a + b + c) and a > c. 
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The strictly positive and convergent invariant measure x is given by 



Xk = T) k (y/ac/a) if X = v, (5) 

= n{[(X-b + VA)/(2a)] fc - [(A -6- VA)/(2a)] fe } ifX>v, (6) 

for k > 1, where n is an arbitrary constant and A = (6 — A) 2 — 4ac. 

Proof. Let M( fe ) denote the fc x fc north-west truncations of M. Then, by a modifi- 
cation of van Doom et al. |15l Theorem 1, Eqn (9)] we obtain 

sp(M^) = min max < b + u i+1 + — \. 

u>0 l<i<h y Ui J 

Then, by Seneta [H Theorem 6.8] 

v = lim sp(Af' fe ^) = min sup \ b + Uj+i + — > = 6 + 2y / ac, 

A:— >oo u>0 j (_ J 

with the last equality following from arguments analogous to those in the proof of 
Theorem 3.2 in Latouche et al. [6]. 

Now, for any x to satisfy xM = Xx, its elements have to satisfy the constraints 

bx\ + 0x2 — Xxi (7) 
axk+i + {b — X)Xk + cxk-i = for k > 2. 

Let A = (b — A) 2 — 4ac. There are three cases, for each of which x takes a specific form 
(Korn and Korn Chapter 20, Section 4.5]). 

Case 1: A = 0. Then, A = b ± 2Vac, and for k > 1, 

x k = a [(A - b)/(2a)] k + c 2 k[(X - b)/(2a)] k , (8) 

where c\ and C2 are constants. Substituting ([8} into ([7]) gives us c\ = 0. To ensure 
Xk > for all k, it is necessary that A > b. Consequently, A = b + 2^/ac, and we obtain 
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For x to be convergent, we require that ^fac/a < 1 and thus a > c. 
Case 2: A > 0. Then, A < b - 2^/ac or A > b + 2 v / ac, and for k > 1, 
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*fc = c 3 [((A - b + VA)/(2a)] k + c 4 [((A - b - VA)/(2a)] k , (9) 

where C3 and C4 are constants. Substituting © into (|7J| gives us C3 = — C4, and ([9]) 
simplies to It is clear from ^ that a; > if and only if A > b. Thus, A > b + 2y/ac. 

Finally, xl < 00 if and only if ((A — b) + \fK)/{2a) < 1, the latter being equivalent 
to A < 2a + b and A < a + b + c. As b + 2^/ac < A < 2a + b. both a < c and a = c lead 
to a contradiction. Consequently, a > c. 

Case 3: A < 0. Then, b ~ 2y/ac < A < b + 2^/ac and 

Xk = (c/a) fe (c 5 cos(fc0) + C6 sin(&0)), (10) 

where = arccos(&/(2-y/ac)), < < tt and C5 and cq are constants. 

Since we are looking for x > 0, Case 3 is not feasible. Indeed, we can rewrite (|TU1) 
as, for fc > 1, Xk = C7(c/a) k cos(cs + kip) where < c$ < 2ir and C7 is arbitrary. It can 
be easily shown that there exists kg such that cos c$ and cos(cs + fco<^) have different 
signs. 

Among the progeny distributions that may be associated with the mean progeny 
matrix given in ([4]), we choose 

Pi{s) = {b/t)s\ + {c/t)s\ +1 + 1 - (6 + c)/t for i = 1, 

= (o/uJaJL! + (&A0 S ™ + (c/u)s^ +1 + 1 - (a + b + c)/u for i > 2, 

where t = \b + c] + 1 and w=[a + 6 + c]+l. By varying a, 6, and c, we shall cover 
the three possible cases q = q<l, q = q = l, and q < q = 1. 
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Figure 1: Case 1. a — b = 1/2, and c = 1/3. Left: the values of q[ k ' (continuous line) and 
of q[ k ^ (dashed line). Right: first entries of q' 20 - 1 (continuous line) and q' 20 ' (dashed line). 

Case 1: q = q < 1. Take a = b = 1/2, and c = 1/3. With these, v = 1.28 > 1, and 
q = q < 1 by Lemma 12.11 and Proposition 14.11 

We illustrate in Figure Q] the convergence of the sequences {q^} and {q^}. On 
the left, we plot q[ k ^ and q[ k ^ for k = 1 to 20; the two sequences rapidly converge to 
a common value approximately equal to 0.89. On the right, we plot qf ^ and qf" ^ 
for i = 1 to 20; we observe that the first 15 entries are well-approximated after 20 
iterations but the next entries require more iterations because for high values of i, the 
approximation process for qi and qi starts with a higher value of k. 

A sequence {xk}k>o converges linearly to x if there exists < /i < 1 such that 
limfe-^oo | a; — a;fc+i|/|a; — Xk\ — fJ>, and /x is called the convergence rate. Our numerical 
investigations indicate that the convergence of q^ as well as that of q^ is linear, for 
fixed i. We give one example in Figure^ where we plot the ratios \qi — q[ k+1,> \/\qi — q[ k ^\ 
and |<7i — q[ k+1 ^ \/\qi — q[ |; not knowing the value of either qi or q\, we have used the 
values at the 20th iteration. The two sequences are seen to converge linearly at the 
same rate \i = 0.26 approximately. 

Case 2: q = q = 1 The parameters here are a = b = 1/2, and c = 1/25. Here a > c 
and for any individual, the type of its descendants drifts over successive generations 
toward type 1, the least prolific of types. The convergence norm is v = 0.78 < 1, which 
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Figure 2: Case 1. a = b = 1/2, and c = 1/3. Convergence rates of the sequences q[ k ^ 
(continuous line) and (dashed line). 

implies that q — 1. We shall conclude from Proposition 15.11 and Proposition 14.51 (with 
A = v) that q = 1 as well, once we show that the dichotomy property holds. The 
progeny generating function is given by 

Pi(s) = (l/4)s? + (l/50)s? +1 + (73/100) for i = 1, (11) 

= (V6)4-i + (1/6)4 + (V75)4 + i + (49/75) for i > 2. (12) 

To verify that the dichotomy property holds, we use the sufficient condition ([3]). In 
view of dm [TJ), we observe that for all i, F[\Z 2 \ = | ip a = i, 1 < \Zi\ < k] > 
(min(73/100, 49/75)) fc , and we conclude that is satisfied with = 2 and = 
(49/75) fc . 

To illustrate the observation made in Subsection 14.41 about the effects of the initial 
type's distribution, we take the parameters 

Ai=^ = 0.78<l and A 2 = 1.02 < 1.04 (= a + b + c). 

By Proposition 15.11 there exist oli > and a.2 > such that a.\M — Xioci and 
ceil = 1, and such that oliM — \2Oi2 and ct2l = 1- 

If ipo has the distribution ai, then lim„E[|Z„|] = 0, while if it has the distribution 
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Figure 3: Case 2: a = b = 1/2, c = 1/25. The first 50 components of the initial type 
distribution vector ati, for which lim n E[j.Z n |] = 0, and of 0:2, for which lim n E[|.2 n |] = 00. 



a.2, then lim n E[|Z„|] = 00. In both cases, extinction is with probability 1. We plot in 
Figure [3] the first 50 components of oli and a.2- The difference between the two is that 
the distribution ol\ is concentrated on small types, so that the process has less chance 
of building a high population before its eventual extinction. 



Case 3: q < q = 1. Take a = 1/25, b = c = 1/2. Here, a < c and v = 0.78 < 1 < 
a + b + c; thus, q < q = 1 but we do not know if q = 1 or not. 

We show on the left of Figure [3] the values of q\ and q\ ' for k = 1 to 60. Judging 
from this, we conclude that q\ < 1 = gj. On the right of that figure, we give <?, and 
qf for i = 1 to 60. 

To confirm the conclusion that q < q = 1, we have simulated the branching process 
and we give one particular sample path on Figure [S] the whole population \Z n \ seems 
to grow without bounds, while individual types appear, grow in importance, and 
eventually disappear from the population. 
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Figure 4: Case 3: a — 1/25, b = c = 1/2. Left: the values of (continuous line) and of 
?i (dashed line). Right: first entries of q' 60 ' (continuous line) and q' 60 ' (dashed line). 



450 
400 
350 
n 300 
| 250 
1200 

Q. 

£ 150 
100 
50 
0, 



type 10 




---type 20 




' - - type 30 




, —type 40 




--type 50 




type 60 


f 


- — total size 






f + /^. : 




• TM h ..... .11- 



20 40 60 80 100 120 140 160 180 
Generation 



Figure 5: Case 3. a = 1/25, b — c — 1/2. Simulation of the evolution of the population size 
in different types and of the total population size 

5.2. Reducible example 

Consider the mean progeny matrix with the structure 



M 



h ci 

hi c 2 

^3 c 3 



where hi > and a > for all i. 
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In this special case of reducible mean progeny matrix we may associate another 
interpretation to the sequence {q^}. Let us define the local extinction of a specific 
type. This event is = [lirrin-joo Z n k = 0], independently of the other types. 

(k) 

A moment of reflection shows that E). = Hi<kEe and, furthermore, that ql is the 
probability that type k eventually becomes extinct, given that the process starts with 
a first individual of type i. This allows us to give another proof that q^ > g( fe+1 ) and 
that the sequence converges to q: 



lim q (k) = lim P[£ fc M = Km P[n^< fe £^ ] 

k— j-oo k—too k— ¥oc 

= V[r\e<<x>Ee ip ] = q- 



In the reducible case, the equation s — P(s) may have more than two distinct solutions 
and, in particular, it is possible that q < q < 1, as we show on one example. 

Take hi — and Cj = 1.9 for every i except for i = 10, where 610 = 1.6 and cxo = 0.8. 
That is, in general, type % ^ 10 individuals have only children of the next type, slightly 
less than two on average, and type 10 is different. If it were not for type 10, the whole 
population would behave as a supercritical process, with each type getting extinct 
after one generation. Individuals of type 10 do reproduce themselves, in a supercritical 
fashion. 

Assume that the progeny generating function is 

P t (s) = (19/30)4fi + (H/30) for z ^ 10, 

= (2/5)4 + (1/5)4+1 + (2/5) for i = 10. 

As the sequence {sp(M^ fe ^)} converges to v — 1.6 > 1, we know by Proposition 14. 21 that 
q ^ 1. Furthermore, Proposition ^. 31 implies that gj < 1 for 1 < i < 10, and qi = 1 for 
i > 11. 

We show {q^} and {q^} on the left in Figure [5] and the plot clearly makes it 
appear that q% < q$ < 1. On the right, we give the values of <7 l - 3 °' ) and q^ 30 ^ for 
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Figure 6: Left: the values of gg fc ' (continuous line) and of (dashed line). Right: first 
entries of q' 30 ' (continuous line) and q*- 30 ' (dashed line). 

1 < { < 30. For i > 11, local extinction has probability 1 since every type exists for 
one generation only, and the global probability, at least if i is sufficiently smaller than 
30, is close to 0.41, the extinction probability of a single-type branching process with 
progeny generating function 



P(s) = (19/30)s 3 + (11/30). 



We thus see that if extinction happens in the single-type process, then it does so in a 
few generations. 
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